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ABSTRACT 

We investigate the apparent discrepancy between gas and dust outer radii derived from millimeter 
observations of protoplanetary disks. Using 230 and 345 GHz continuum and CO J=3-2 data from 
the Submillimeter Array for four nearby disk systems (HD 163296, TW Hydrae, GM Aurigae, and 
MWC 480), we examine models of circumstellar disk structure and the effects of their treatment of 
the outer disk edge. We show that for these disks, models described by power laws in surface density 
and temperature that are truncated at an outer radius are incapable of reproducing both the gas and 
dust emission simultaneously: the outer radius derived from the dust continuum emission is always 
significantly smaller than the extent of the molecular gas disk traced by CO emission. However, a 
simple model motivated by similarity solutions of the time evolution of accretion disks that includes a 
tapered exponential edge in the surface density distribution (and the same number of free parameters) 
does much better at reproducing both the gas and dust emission. While this analysis does not rule 
out the disparate radii implied by the truncated power-law models, a realistic alternative disk model, 
grounded in the physics of accretion, provides a consistent picture for the extent of both the gas and 
dust. 

Subject headings: accretion, accretion disks — circumstellar matter — planetary systems: protoplan- 
etary disks — stars: pre-main sequence 



1. INTRODUCTION 

Characterizing the gas and dust distribution in the 
disks around young stars is important for understand- 
ing the planet formation process, as these disks provide 
the reservoirs of raw material for nascent planetary sys- 
tems. A common method of modeling circumstellar disk 
structure is to use models described by power laws in 
surface density and temperature that are truncated at 
a particular outer radius. This prescription has its his- 
torical roots in calculations of the minimum mass so- 
lar nebula, which indicated a surface density profile of 
E oc r~ 3 / 2 (e.g. IWeidenschillinel fl977l) . as well as the- 
oretical predictions of a radial power-law dependence 
of temperature for accreting disks arou nd young stars 
([Adams fc Shvj||1986l lAdams et alJll987[ ). Observation- 
ally, the parameterization of temperature and surface 
density as power-law functions of radius began with early 
spatia lly unresolved studies of continuum emission from 
disks (jBeckwith et all 119901: iBeckwith fc Sargenti 119911 ). 
These models have since been refined and applied to spa- 
tially n3sj3Reji_j3bser\^^ 
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1998), and they have proven useful for 



understanding the basic global properties of disk struc- 
ture. Recently, however, with the advent of high signal- 
to-noise, multi-frequency observations of gas and dust 
in protoplanetary disks, these models have begun to en- 
counter difficulties, particularly in the treatment of the 
outer disk edge. 

The extent of the gas and dust distribution in circum- 
stellar disks has implications for our understanding of the 
planet formation process in our own solar system. There 
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is some evidence for a sharp decrease in the surface den- 
sity of Kuiper Belt objects beyond a distance of 50 AU 
from the Sun dJewitt et alj |1998t iTruiillo fc BrowiJl200H 
iPetit et ail l2006h. However , the origin of this edge is 
unclear. lAdams et al.l (|2004f ) note that the observed dis- 
tance is far interior to the radius at which truncation 
by photoevapo r ation would be expected to occur, while 
lYoudin fc Shi (120021 ) find that the presence of such an 
edge in planetesimal density could be explained by drift- 
induced enhancement. A compelling possibility is that 
the Sun formed in a cluster environment, and the early 
solar disk wa s truncated by a close encounter with a pass- 
ing star (see lReipurth| [2005. and references therein). A 
more complete understanding of the outer regions of pro- 
toplanetary disks may provide insight into the processes 
t hat shape th e oute r solar system. 

iPietu et al.l (|2005f ) present multiwavelength millimeter 
continuum and CO isotopologue observations of the disk 
around the Herbig Ae star AB Aurigae and found from 
fitting models of disk structure described by truncated 
power laws that the outer radius of the dust derived 
from continuum emission (350±30 AU) was much smaller 
than that of the gas derived from 12 CO J=2-l emission 
(1050± 10 AU). They suggest that a change in dust grain 
properties resulting in a drop in opacity could be respon- 
sible for the difference, and note the possible association 
with a ring featu re in the disk at 20 AU. A similar result 
was obtained bv llsella et al.l (|2007f ) from observations of 
the disk around the Herbig Ae star HD 163296: they 
found a significant discrepancy between the outer radius 
derived for the dust continuum emission (200 ± 15 AU) 
and that derived from CO emission (540 ±40 AU). These 
data appeared to require a sharp drop in surface density, 
opacity, or dust-to-gas ratio beyond 200 AU; however, as 
they discuss, there is no obvious phys ical basis for such 
a discontinuity. As llsella et al.l (120071 ) demonstrate, the 
discrepancy in outer radii derived from the dust and gas 
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is not simply an issue of sensitivity; the observations were 
sufficiently sensitive to detect emission from the power- 
law dust disk if it did extend to the radius indicated by 
the CO emission. The underlying issue is that the trun- 
cated power law model does not simultaneously repro- 
duce the extent of both the continuum and CO emission 
for these disks. 

Using data from the Submillimctcr Array 3 we show 
that the same apparent discrepancy in gas and dust outer 
radius applies to the circumstellar disks around several 
more young stars. In an attempt to understand the origin 
of this discrepancy, we investiga te an alternative surfac e 
density profile based on work by Hartma nn et al.l (|1998f) , 
which is similar to a power law profile in the inner disk 
but includes a tapered outer edge. We show that this 
model, which has a physical basis in similarity solutions 
of disk evolution with time, is capable of simultaneously 
reproducing both continuum and CO emission from these 
disks. The primary difference between this model and the 
truncated power-law disk is that instead of a sharp outer 
edge the surface density falls off gradually, with sufficient 
column density at large radii that CO emission extends 
beyond the point at which dust continuum emission be- 
comes negligible. 

2. MILLIMETER/SUBMILLIMETER DUST CONTINUUM 
AND CO J=3-2 DATA 

The analysis was conducted on extant SMA data of 
the disks around of HD 163296, TW Hydrae, CM Au- 
rigae, and MWC 480. The dates, frequencies, antenna 
configurations, number of antennas, and original publi- 
cations associated with the data sets are listed in Table 
Q] The four disk systems chosen for this analysis are all 
nearby, bright, isolated, and have been well studied at a 
wide range of wavelengths. The velocity fields of these 
disks all appear to be well described by Keplerian rota- 
tion (llsella et al JI2007I : lOi et al.|[200l iDutrev et~aLlll998t 
iPietu et al.ll2007t ). The relevant properties of these sys- 
tems (spectral type, distance, stellar mass, age, and disk 
inclination and position angle) are listed in Table [2] 

3. DISK MODELS 

Using the SMA data available for the four disk systems, 
we compared two classes of disk models: the first model 
is described by power laws in surface density and temper- 
ature and is truncated at an outer radius R ou t (details 
in §3.1|) . and the second model is described by a power 
law in temperature and a surface density profile similar 
to a power law in the inner disk but tapered with an ex- 
ponential edge in the outer disk (details in ^3.2p . This 
latter model is not intended to be a definitive description 
of these disks, but rather illustrative of the broader cat- 
egory of models without a sharp outer edge. The model 
fitting process involved deriving a minimum \ 2 solution 
for those parameters of each class of model that best 
fit the continuum emission, and then using standard as- 
sumptions to predict CO emission (described in i)3.4p . 
The CO emission was not used to determine the model 
fits, due to the computational intensity of solving the ex- 
citation and radiative transfer for the molecular line for 

3 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



a large grid of models. 

3.1. Truncated Power Law 

For the t runcated power law models, we used the pre- 
scription of IDutrev et al.l (|1994h . In this framework, the 
disk structure is described by power laws in tempera- 
ture and surface density, with the scale height specified 
through the assumption that the disk is in hydrostatic 
equilibrium: 



T(R) = T wo 
E(i?) = E 100 
H(R) = 4/ 



R 



1QQAU, 

f R x ~" 
\100AU 

l 2R*k B T k {R) 



(1) 
(2) 
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where the subscript TOO' refers to the value at 100 AU, 
ks is Boltzmann's constant, G is the gravitational con- 
stant, M+ is the stellar mass, and mo is the mass per 
particle (we assume 2.37 times the mass of the hydrogen 
atom). Combining these expressions and the assumption 
of hydrostatic equilibrium, the volume density n(R, z) is 
given by: 



n(R, z) = 



g(g) 

y/*H(R) 



cxp -{z/H(R)f 



(4) 



where z i s the vertical height above the midplane. As 
noted bv IDutrev et al.l (|2007l ). this definition implies a 
scale height of H(r) = v^Cs/fi, where c s is the sound 
speed and Q the angular velocity, while other groups use 
H(r) — c s /0; this difference should be taken into ac- 
count when comparing our results with other disk struc- 
ture models. During the modeling process, we recast the 
surface density normalization in terms of the midplane 
density at 100 AU, so that the parameter Eioo is replaced 
by fT-ioo- This power-law model of disk structure has five 
free parameters: Tioo, 9, "ioOj P, and Rout- 

3.2. Similarity Solution from Accretion Disk Evolution 

While versatile and ubiquitous, the truncated power 
law models of disk structure have one obviously unphys- 
ical feature: a sharp outer edge. In the absence of dy- 
namical effects (e.g. from a binary companion) or large 
pressure gradients to confine the material, disk struc- 
ture at the outer edge is expected to taper off grad- 
ually. A description of the structure of an isolated, 
ste adily accreting di s k as i t evolves with time is provided 
bvlHartmann et al.l (Il998f). w ho expand on the work of 
iLvnden-Bell fc Pringld (|1974f ) to show that if the viscos- 
ity in a disk can be written as a time-independent power 
law of the form v cx i? 7 , then the similarity solution for 
the disk surface density is given by 



E(r) = — T 



-(5/2-7)/(2-7) 



exp 



,2-7 



(5) 



where C is a constant, r is the disk radius in units of 
the radial scale factor i?i such that r = R/Ri, and T 
is the nondimensional time T = t/t s + 1 where t ? is th e 
viscous scaling time (eq. 20 in lHartmann et ail [l998V 
For simplicity, when applying these models to our data 
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TABLE 1 

Sources of SMA 230/345 GHz continuum and CO J=3-2 data. 



Freq./ 








No of 




TV (i i is it i o n 


Object 


D t 


Config. 


Antennas 


Reference 


230 GHz 


HD 163296 


23/24 Aug 2003 


Compact N 


7 


i 




TW Hydrae 


10 Apr 2005 


Extended 


8 


2 






27 Feb 2005 


Compact 


8 


2 




GM Aurigae 


10 Dec 2006 


Extended 


8 


3 




MWC 480 


18/20 Nov 2003 


Compact N 


8 


1 


345 GHz/ 


HD 163296 


23 Aug 2005 


Compact 


8 


4 


CO J=3-2 


TW Hydrae 


28 Dec 2006 


Compact N 


8 


5 




GM Aurigae 


26 Nov 2005 


Compact 


7 


6 




MWC 480 


21 Oct 2005 


Compact 


8 


1 



References.— (1) SMA a rchive; (21 101 et all pOOlih ; (3) 
(6) lAndrews fc Williams! J2007h . 



l'i et al. (in prep); (4) Isclla ct al. (2007) ; (5) Qi et al. (2007, submitted); 



we used physical units and absorbed several of the pa- 
rameters into two constants so that the surface density 
is of the form 



£(i?) = -^exp 



2-7 



(G) 



where R is the disk radius in AU and ci, C2, and 7 are 
constants that we allowed to vary during the fitting pro- 
cess. 

The temperature profile for the similarity solution disk 
model is identical to that of the truncated power-law 
disk, except that its spatial extent is infinite. We do not 
allow it to drop below 10 K, but this low temperature 
limit is not problematic for any of the disks considered 
here. This model therefore includes five free parameters: 
Tioo, q, ci, 7, and C2. The constant c\ describes the 
normalization of the surface density similar to n.100 in 
the power-law model, while the constant C2 is analogous 
to the outer radius, since it describes the radial scale 
length over which the exponential taper acts to cause 
the surface density to drop towards zero. 

3.3. Model Comparison 

The surface density description for the similarity solu- 
tion is similar to the truncated power law except at the 
outer edge of the disk. In the inner regions of the disk for 
which i? < C2, we may expand the exponential so that 



exp [~(R/c 2 ) 2 -'\ 
density becomes 

E(J2) = -^(1 



1 - 

V Co. ' 



_R\2- 7 



2-7 



SL _ _^R2(l-7) (7) 



and the surface 



Cl 



In the a- viscosity context (|Shakura &; Svunvaevlll973T) . 
for a vertically isothermal disk with the typical temper- 
ature index q = 0.5, we would expect that 7 = 1. This 
implies that for standard assumptions, the inner disk sur- 
face density will be described by a power law in R with 
index 7, modified by a constant (^) due to the influence 
of the exponential. If 7 deviates from 1, an additional 
shallow dependence on R would be expected. 

It is illuminating to consider the behavior of these mod- 
els in the Fourier domain, the natural space for interfer- 
ometer observations. To do so we define the coordinate 
1Z UV , the distance from the phase center of the disk in 
the (u, v) plane, as it would be observed if the disk were 
viewed directly face-on. To perform the deprojection 



from the inclined and rotated sky coordinates, we cal- 
culate the position of each point in the (it, v) plane as a 
projected distance from the major and minor axes of the 
disk, respectively: d a = 72sin</> and db = 72cos</>cosi, 
where i is the disk inclination, 72 = (u 2 + v 2 ) 1 / 2 , <fi 
is the polar angle from the major axis of the disk, 
(p = arctan(?;/u— PA), and PA is the position angle mea- 
sured east_ofjior^_yien 1Z UV = (d 2 + d 2 ) 1 / 2 (|Lav et al.l 
ll997tlHughes et al.ll200l . 

In the Fourier domain, the truncated disk models show 
"ringing" and the visibilities will drop rapidly to zero in 
the vicinity of 1Z UV = l/R ou t- Under the simplifying 
assumption of 7 = 1, the Fourier transform of the sim- 
ilarity solution becomes a convolution of two functions: 
(1) 1/72., which is just the Fourier transform of the 1/R 
dependence of a power law extending from zero to infin- 
ity, and (2) c 2 /(l + 72 2 c 2 ) 3 / 2 , where c is a scaling con- 
stant for the term that describes the exponential taper. 
Since these two functions both decrease monotonically 
and are always positive, the visibilities drop smoothly to 
zero without any ringing. 

3.4. Model Fitting 

For both disk models, we fit for the five parame- 
ters describing the disk temperature and surface den- 
sity structure using the continuum data for each disk 
with the widest range of available baseline lengths. The 
position angle and inclination were fixed and adopted 
from previous studies (see Table [2]). For opacity, 
we assume the stan d ard m illimeter value adopted by 
iBeckwith fc Sargentl (| 19911 ) (k v = k (v/vq) 13 , where 
«o = 10.0 cm 2 /gd U st, Wo=l THz, and (3 — 1), although 
we allow [3 to vary in order to obtain the proper normal- 
ization when extrapolating from one frequency to an- 
other. Due to the ~ 100 AU spatial resolution, these 
data are not sensitive to the inner radius of the disk. For 
both sets of models, therefore, we sim ply fix the inner ra- 
dius at a value of 4 AU for TW Hya jCalvet et al1l2002t 
Hughes et al]l2007h . 24 AU for GM Aur (jCalvet et all 



20051 ) and 3 AU for the other two systems, for which 
reliable inner radius information is not available; this 
is sufficiently small that changes in the inner radius do 
not affect the derived model parameters. To compare 
the models to the data, we use the Monte-Carlo radia- 
tive t ransfer code RATRAN ( H ogerheiide fc van der Takl 
2000) to calculate sky-projected images of the dust con- 
tinuum and CO emission, with frequency and bandwidth 
appropriate for the observations, and assuming Keple- 



4 



rian rotation. We then use the MIRIAD task uvmodel to 
simulate the SMA observations, with the same antenna 
positions and visibility weights. 

For each set of parameters, we directly compare the 
model visibilities to the continuum data and calculate 
a x 2 value, using the minimum x 2 value to determine 
the best-fit parameters. The resulting best-fit models 
are shown along with continuum data for both frequen- 
cies in the left panels of Figure [TJ The abscissa gives 
the deprojected radial distance in the (it, v) plane, and 
the ordinate shows the real and imaginary components 
of the visibility. For a disk with circular symmetry, the 
imaginary components should average to zero. The 230 
GHz continuum data are depicted as open circles, while 
the 345 GHz data are filled circles. The best-fit power- 
law model is shown in blue and the similarity solution in 
orange. Dotted and dashed lines distinguish between the 
230 and 345 GHz model predictions; the fit was deter- 
mined at that frequency with the largest baseline cover- 
age and extrapolated to the other frequency by varying 
0. The uncertainties quoted for /3 reflect an assumed 
10% calibration uncertainty. Note that varying (3 has no 
effect on the modeled CO emission. 

We measure values of f3 that are consistent with 1, 
which is in agreement with the typical values mea- 
sured for disks in the Taurus-Auriga as sociation (e.g. 
iDutrev et all 119961 : iRodmann eFail I2006D . These shal- 
low millimeter spectral slopes indicate that some grain 
growth has occurred from ISM grain sizes, which typi- 
cally exhibit a steeper spectral slope (f3 ~ 2). In partic- 
ular, the value of 1.2 measu red for GM Aur matches we ll 
the value of 1.2 reported by I Andrews fc Williams! {2007), 
and the value of 0.7 measured for TW Hydrae ma tches 
well the value of .7 reported bv lCalvet et alj (|2002| ) and 
iNatta et al.l (|2004h . 

After fitting the continuum, we then assumed a 
gas/dust mass ratio of 100 and a standard interstel- 
lar CO/H2 mass ratio of 10 -4 to predict the expected 
strength and spatial extent of CO emission from the 
disks, based on the best-fit continuum model. We as- 
sume throughout that the gas and dust are well-mixed, 
and that CO traces molecular hydrogen. We do not take 
into account the complexities of disk chemistry, such as 
the de pletion of CO molecules in the co ld, dense mid- 
plane |Aikawall20"07l ISemenov et alj |2006). However, de- 
viations from these simple assumptions should have no 
appreciable effect on our conclusions concerning the ra- 
dial extent of CO emission. 

Since we neglect the vertical temperature gradient in 
the disk, we might expect to underpredict the strength 
of optically thick CO line emission, which likely origi- 
nates in the upper layers of the disk that are subject to 
heating by stellar irradiation. The continuum emission, 
by contrast, is likely weighted toward the cold midplane 
of the disk. For this reason, after obtaining an initial fit 
from the continuum, we allowed the temperature scale 
(Tioo) to vary to best reproduce the flux levels of the ob- 
served CO emission, and then iteratively fit for the other 
structural parameters (q, n\oo, Cz/Routi and j/p)- 

Deriving the temperature from the CO emission in this 
way may underestimate the midplane density in some 
cases, due to the degeneracy between Tioo and nioo: 
the temperature derived from CO emission is typically 
greater than or equal to that of the shielded midplane, 



depending in detail on the dust opacity and molecular 
dissociation due to ultraviolet radiation in the upper 
disk layers (for a discussion of the proce sses involved, see 
I Jonkheid etail l2007t llsella et al.l I2007D . For the disks 
considered, the temperature derived for the dust contin- 
uum emission was within ~ 40% of that derived to match 
the CO line strength. 

4. RESULTS AND DISCUSSION 

The parameters for the best-fit model solutions to the 
continuum data for each source and for each of the two 
model types are listed in Table [3j This table lists only 
the set of parameters with the minimum x 2 value; for- 
mal errors are not quoted as these are not intended to 
be definitive structural models but simply illustrative of 
the differences between the model classes in their treat- 
ment of the outer edge. The midplane surface density 
profiles for these models are plotted in Figure [2J The 
solid lines depict the profile for the power-law solution, 
while the dashed lines are for the similarity solution. The 
parameters of the two model solutions are very similar, 
particularly for HD 163296 and MWC 480. For all four 
disks, the two model solutions are particularly similar 
just within the outer edge of the disk, around the range 
of radii well-matched to the resolution of the data (~ 200 
AU for HD 163296, ~ 90 AU for TW Hydrae, ~ 200 AU 
for GM Aurigae, and - 300 AU for MWC 480). The 
outer radius for the power-law solution typically falls at 
roughly twice the scale length (02) of the similarity solu- 
tion. The analogous parameters 7 and p, which describe 
how quickly the midplane density drops with radius, are 
also very similar between the two models. 

CO J=3-2 emission predicted from these best-fit mod- 
els is shown in the right panel of Figure [TJ The similar- 
ity solution is shown in the blue-contoured central plot, 
and the power-law model in the orange-contoured plot 
on the right. Recessional velocity is plotted on the ab- 
scissa while the position offset along a slice through the 
disk major axis is plotted on the ordinate. The horizon- 
tal dashed line in each figure represents the extent of the 
outer radius (R ou t) derived for that source in the context 
of the truncated power-law model. For all four sources, 
the extent of molecular gas emission from the similarity 
solution is much more closely matched to the data than 
that of the power-law model, even though both reproduce 
the continuum dust emission equally well. 

From Figure [lj it is clear by eye that for all four 
sources, the extent of the CO emission is severely under- 
predicted by the power law model but matches well the 
predicted emission from the similarity solution model. 
A calculation of the x 2 value comparing the predicted 
CO emission for the two models to the observed emis- 
sion shows that the similarity solution matches the data 
better than the truncated power-law model for all of the 
disks in our study. The difference is at the 2a level for 
MWC 480, for which there is only short-baseline data 
with relatively low signal-to-noise, and at the 4tr level for 
GM Aur; for TW Hydrae and HD 163296, the x 2 analy- 
sis shows that, formally, the similarity solution provides 
a better fit to the CO emission than the power-law model 
at the > lOer level. 

The tapered edge of the similarity solution density dis- 
tribution evidently permits a large enough column den- 
sity to produce detectable CO 3-2 line emission, even 
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TABLE 2 
Stellar and disk properties 





Spectral 


Dist. 


Stellar 


Age 


Disk 


Disk 


System 


Type 


(PC) 


Mass (M Q ) 


(Myr) 


PA (°) 


in 


HD 163296 


A1V 


122 a 


2.3 a 


5 fc 


145 c 


46 c 


TW Hydrae 


K8V 


51 d,e 


0.6 


5-20^9 


-45 h 


7'' 


GM Aurigae 


K5V 


140 


0.8* 


2-10^ fe 


51* 


56*' 


MWC 480 


A3 


140 m 


1.8" 


7-8 n '° 


143P 


37? 



Ref eren ces. — (a) Ivan den A nckcr ct al. (1998); (6) Natta ct al. (20 QD; (e) IIsella et al] (t20O7T); (d) Mamaick (2005); (e) Hof f et al.l 
(Il998l); (f) IKastner et al.l (119971); (g) I Webb et all il99SlV; (h ) IQi et al.l 1120041); (i) IDutrev et al.| 111998); (7) IBeckwith et all 11990) ; ffc) 
ISimon fc Pratol I1995fi 77ro) Ifhe et aLI (119941) : fn) IPietu et all (120071) ; (o) ISimon et all l|200a) : M iHamidouche et al.l 1 120061) . 
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1.0 
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similarity, power low 



230 GHz 
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Fig. 1. — Comparison between the data and the two types of models (similarity solution and power law) for the four disks in our 
sample: (a) HD 163296, (b) TW Hydrae, (c) GM Aurigae, and (d) MWC 480. For each source, the left panel shows the real and imaginary 
visibilities as a function of deprojected (if, v) distance from the phase center. Symbols are SMA data; open circles are 230 GHz and filled 
circles are 345 GHz continuum. The lines represent the best fit to the 345 GHz continuum for the power law (orange) and similarity (blue) 
models. Dashed lines show the model at 345 GHz while solid lines are 230 GHz. The right panel shows position-velocity diagrams of the 
J=3-2 rotational transition of CO along the major axis of the disk. The left plot (black contours) shows the SMA data. The middle plot 
(blue contours) displays the emission predicted by the similarity solution parameters that provide the best fit to the continuum emission, 
and the right plot (orange contours) displays the emission predicted for the best-fit power-law model. The horizontal dashed line across 
the right panel represents the extent of the outer radius (Rout) derived for each source through fitting of the continuum emission in the 
context of the power-law model. The contour levels, beam, and velocity resolution for each source are as follows: (a) [2,4,6,8,10,12] X 1.1 
Jy/beam, 3.0x2.1 arcsec at a position angle of 14.3°, and 0.35 km/s; (b) [2,4,6,8] x2.0 Jy/beam, 4.0x1.8 arcsec at a position angle of 3.2°, 
and 0.18 km/s; (c) [2,4,8,12,16] x0.5 Jy/beam, 2.3x2.1 arcsec at a position angle of 12.9°, and 0.35 km/s; (d) [2,4,6,8,10] x0.5 Jy/beam, 
2.5x2.3 arcsec at a position angle of 45.3°, and 0.35 km/s. 



though it has dropped off enough that the continuum 
emission is negligible. The power-law model, by con- 
trast, is strictly limited in the extent of its CO emission 
by the sharp outer radius. In particular, for the case of 
HD 163296, the CO emission predicted by the power law 
model (orange contours in the right panel of Figure [T^,) 
falls to 4.4 Jy/beam at a distance of 1.8 arcsec (220 AU) 
from the source center, while the similarity solution (blue 
contours) maintains this brightness out to a distance of 



4.7 arcsec (570 AU). This latter size is well matched to 
the data (black contours) which extends at this bright- 
ness to a distance of 5.0 arcsec (600 AU). These dis- 
tances likely overestimate the true physical extent of the 
disk due to convolution with the 2.1 x 3.0 arcsec beam, 
th ough they are very comparable to the values observed 
bv llsella et al.l (|2007h . While the similarity solution does 
not provide a perfect fit to the data, nor do we intend 
it to do so, it illustrates that the outer radius discrep- 
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TABLE 3 

Parameters for. best-fit continuum models 



Source 


Model 


x 2 


Tioo (K) 


g 


ni a (cm -3 ) 


c 2 (AU) 
Rout (AU) 


7 
P 





HD 163296 


Similarity 
Power Law 


2.29 
2.26 


65 
60 


0.4 
0.5 


5.3 x 10 11 
6.7 x 10 11 


125 
250 


0.9 
1.0 


u - 4 -0.3 

5 +0 ' 5 
u -°-0.3 


TW Hydrae 


Similarity 
Power Law 


2.42 
2.41 


40 
30 


0.2 
0.5 


2.3 x 10 11 
7.1 x 10 10 


30 
60 


0.7 
1.0 


u - '-0.1 
u - '-0.1 


GM Aurigae 


Similarity 
Power Law 


2.19 
2.17 


50 
40 


0.5 
0.4 


1.1 x 10 11 
5.0 x 10 11 


140 
275 


0.9 
1.3 


i o+U.b 
± * z -0.1 
1 q+0.5 


MWC 480 


Similarity 
Power Law 


1.86 
1.86 


50 
45 


0.8 
0.7 


1.0 x 10 11 
1.3 x 10 11 


200 
275 


1.1 
1.3 


-0.4 
n 7+0-5 
u - '-0.4 



ancy is peculiar to the truncated power-law model; other 
disk structure models with a tapered outer edge may be 
able to reproduce the gas and dust emission as well as, 
or better than, the similarity solution adopted here. 

Analysis of the CO excitation in the similarity solu- 
tion model shows that the extent of the CO J=3-2 emis- 
sion in these disks coincides roughly with the radius at 
which the line excitation becomes subthermal, deter- 
mined primarily by where the mid-plane density drops 
below the critical density (~ 4.4 x 10 4 cm~ 3 at 20 K, 
though effectively lowered when photon trapping plays a 
role). In the similarity solution model, the surface den- 
sity distribution steepens dramatically at large radii, but 
without the sharp truncation of the power-law model. 
This suggests that caution should be exercised not only 



when comparing outer radius measurements based on 
dust continuum and molecular gas emission, but also 
when comparing measurements based on emission from 
different transitions of CO or from isotopologues of the 
CO mol ecule that have dif fering abundances and optical 
depths. iPietu et ail (|2007l ) fit truncated power law mod- 
els to the disks around DM Tau, LkCa 15, and MWC 480 
in several different isotopologues and rotational transi- 
tions of CO. For the two cases in which multiple tran- 
sitions of the 13 CO molecule were observed, the derived 
outer radius is marginally smaller for the J=2-l transi- 
tion than the J=l-0 transition. This result is consistent 
with the expected trend that lower-J transitions will ex- 
hibit larger outer radii due to their lower critical density: 
a lower critical density will be reached at a greater dis- 
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Fig. 2. — Midplane density structure of the models that pro- 
vide the best fit to the continuum data. Solid lines show truncated 
power-law models while dashed lines show similarity solution mod- 
els. 

tance as the surface density t apers off nea r the o uter edge 
of the disk. In all cases the iPietu et al.l (|2007l ) analysis 
also yielded a smaller outer radius in 13 CO than in 12 CO, 
as well as a flatter surface density power law index for 
13 CO than for 12 CO. These differences may be related to 
selective photodissociation, or other chemical processes. 
However, the trends of smaller outer radius and shal- 
lower surface density index in 13 CO are also consistent 
with surface density falling off rapidly at large radii, as 
expected for a disk with a tapered outer edge. In the 
similarity solution model, the less abundant 13 CO iso- 
topologue will become undetectable at smaller radii than 
12 CO, which is more sensitive to the exponential drop in 
surface density in the outer disk. 

It is noteworthy that studies of six largest "proplyds" 
with the most distinct silhouettes in the Orion Neb- 
ula Cluster reveal radial profiles in extinction that are 
well-described by an exponenti al taper at the outer edge 
(jMcCaughrean fc O'Dell I1996D . These isolated disks 
may be analogous to the systems considered here with 
a tapered outer edge. 



Models with tapered outer edges also aid in addressing 
discrepancies between the size of the dust disk observed 
in the millimeter and the extent of scattered light ob- 
served in the optical and near-infrared. For example, 
coronographic observations of TW Hydrae detect scat- 
tered light to a distance of ~ 200 AU from the star 
dKrist et alJ 120001 : iTrilling et al.|[200lt IWeinberger et all 
2002), while the truncated power-law model places the 
outer edge of the dust disk c loser to 60 AU. Si milarly, 
observations of HD 163296 bv lGradv et~aTl (|2000D detect 
scattered light out to ~ 400 AU from the star, much 
larger than the 250 AU radius of the dust disk implied 
by the truncated power-law model. While the exponen- 
tial taper causes the density of the similarity solution to 
drop rapidly with radius, these models retain a substan- 
tial vertical column density for several exponential scale 
lengths. It is therefore plausible that scattered light can 
remain visible at this distance, in contrast to the case of 
the smaller truncated power-law disk. 

Although we intend for the similarity solution applied 
here to be an illustrative rather than definitive descrip- 
tion of the disk structure, it is important to note that 
the particular form applied here has potential implica- 
tions for the study of the evolutionary status of these 
dis ks. The form of the si milar ity so lution developed 
by iLynden-Bell k. Pringld (|1974h and lHartmann et all 
(|1998f ) relates the observed structure to the disk age, vis- 
cosity, and initial radius. Although all three of these vari- 
ables are poorly constrained by current observations, a 
large and homogeneous sample of objects studied in this 
way might reveal evolutionary trends in the disk struc- 
ture. 

5. SUMMARY AND CONCLUSIONS 

With the advent of high signal-to-noise interferometer 
observations that resolve the outer regions of nearby pro- 
toplanetary disks, an apparent discrepancy has emerged 
between the extent of the dust continuum and m olecular 
gas emission (jPietu et al.ll2005l: llsella et al.ll2007[ ). Using 
multi-frequency interferometric data from the Submil- 
limeter Array, we have investigated this disparity for four 
disk systems (HD 163296, TW Hydrae, CM Aurigae, and 
MWC 480) in the context of two distinct classes of disk 
structure models: (1) a truncated power law, and (2) a 
similarity solution for the time evolution of an accretion 
disk. The primary difference between these models is 
in their treatment of the disk outer edge: the abruptly 
truncated outer edge of the power-law disk causes the 
visibilities to drop rapidly to zero, leading to an inferred 
outer radius that is small in comparison with the ob- 
served molecular gas emission. The similarity solution, 
by contrast, tapers off smoothly, creating a broader vis- 
ibility function and allowing molecular gas emission to 
persist at radii well beyond the region in the disk where 
continuum falls below the detection threshold. The outer 
radius discrepancy appears to exist only in the context 
of the power-law models. 

In light of this result, it appears that an abrupt change 
in dust properties for these disks is unlikely, as there is 
no physical mechanism to explain such a discontinuity. 
This may imply that a sharp change in dust properties in 
the early solar nebula is similarly an unlikely explanatio n 
for the Kuiper belt edge observed bv lJewitt et alJ (|1998h . 
and that a dynamical mechanism such as truncation by 



a close encounter with a cluster member (Rcipurth l2005l . 
and references therein) may provide a more plausible ori- 
gin. In this case, we would expect to observe disks with 
sharp outer edges only in clustered environments, and a 
model with a tapered edge would be a more realistic pre- 
scription for investigating the structure of a typical iso- 
lated disk. The tapered disk models provide a natural ex- 
planation for the disparate outer radii observed using dif- 
ferent probes of the disk extent, including co mparison of 
continuum and molecu lar gas observations (|Pietu et al.l 
I2005t llsella et al.ll2007lh and also comparison of different 
isotopolo gues and rotation al transitions of a particular 
molecule (jPietu et al.ll2007h . When predicting CO emis- 
sion, this simple model does neglect potential variance 
in the CO abundance due to depletion in the midplane 
and photodissociation at the disk surface; however, the 
results presented are intended simply to illustrate the 
global differences between gas and dust emission from 



the two model classes, independent of detailed CO chem- 
istry. 

While we cannot rule out disparate gas and dust radii 
in these disks, we show that an alternative disk struc- 
ture model, grounded in the physics of accretion, resolves 
the apparent size discrepancy without the need to invoke 
dramatic changes in dust opacity, dust density, or dust- 
to-gas ratio in the outer disk. 
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